UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 

LIBRARY ROUTIN E M l6 - I87 

TITLE Matrix Multiplication When the Common Dimension 

Is Large 

TYPE Complete program 

2 
DURATION n (2^d + .2p) + n (31. 5p + 25d) milliseconds where 

p = the common dimension of the matrices 

n = the number of rows in the premultiplier 

plus the number of columns in the post- 

multiplier 
d = the number of characters in each element 

of the product matrix 

DESCRIPTION Forms the matrix product .C. = .A B. and then 

d k jppk 

punches the matrix .C, by columns, an N following 
the last element of each column. 

CAPACITY j + k < ij-O; no limit on p 

INSTRUCTIONS TO ILLIAC OPERATOR 1. Program tape 

2 . Parameter tape Black up 

3* Data tape Black up 

k. Continue input of program tape Black up 

Ends on OF order from 018. 

PREPARATION OF TAPES The following parameter tape is punched 

where £, 11, .d, and £ are as defined, above. 

Before the elements of the matrices A and 
B are punched on the data tape they must be scaled 
by a constant which will make (a ) all of the 
elements less than unity and (b) each of the 
row sums of squares of the premultiplier A and 
the column sums of squares of the postmultiplier 
B less than ten. (The. program scales each of the 
elements of the product matrix by 10 while 
forming the product.) The scaled elements are 
arranged in the form of an n by £ matrix composed 



EXAMPLE 



-2- 



of the scaled premultiplier followed by the 



A 

« • • » 



B 



.T 



transpose of the scaled postmultiplier 

This composite matrix is punched by 

columns with an N after the last element in each 

of the columns . 

Suppose we want to find the product matrix where 

.6 
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•5 
1.0 

•3 



• 7 

.2 



2\ 



k 



In this case we have to scale the elements of A 

and B by 10*" . The product matrix will be scaled 

-2 
by the square of this, 10 , and an additional 

10 which is added by the code. Since there will 

be three scaling zeros, to get answers to three 

decisel places we will have to take d equal to 6. 

The parameter tape is punched as follows: 
spaces 2S7N6F3L spaces. 

The data tape is punched as follows: 
+05 +10 +03 +7 +3 +1 +6 N 
+06 +07 +02 ++ +5 +2 +9 ir 

The printed results will be: 
+000590 
+000980 
+000290 
N 

+000^50 
+000650 
+000190 

N 



+000170 
+0002^0 
+000070 
N 
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+000840 
+001230 

+000360 

N 
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Routine X 1 - 18 

00 115K I 
; L5 28F j 
I k2 1L l from 3 

j hi 8F ; 
j kl ( )F ; by 
j F5 1L 

k2 1L 

LO 111F 

32 1L 

22 35F 

00 F 

00 10K 

L5 2F from 13 

LO 5F 

40 120F 
26 19L 
00 F 
00 F 
00 F 
00 F 
kl 3F 

41 F 

81 kF from 9 

LO 112F 
32 9L 
hk 112F 
50 F 
jk 112F 
S5 F 
kO F 
26 5L 

42 10L from 6 
L5 F 
kO ( )F by 9 



Decimal Order Input 



— Clear numerical storage locations 



— Form j 



L 



Input parameters 



J 
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F5 3F 
40 3F 






12 


LO 113F 
32 kL 






13 


92 191F 








26 OL 






lk 


00 F 


—i 




00 F 






15 


r\/-\ / Yrn 


by 21 




16 


00 ( )F 
00 ( )F 


by 21 
by 22 


— Dependent parameters storage 




00 ( )F" 


by 22 




IT 


80 F 








00 ( )F 


by 23 . 




18 


00 121F 
00 121F 




Beginning of numerical storage 


19 


50 2F 
00 59F 


from 1 




20 


L4 2F 
kO F 






21 


Lk 18L 








kO 15L 




— Form dependent parameters 


22" 


Lk F 
ko 16L 






23 


L5 2F 
k2 17L 






24 


L5 IF 
kO 9F 






25 


2k 999F 
41 3F 




Stop to read in second part of master tape 


26 


40 121F 
50 26L 






27 


26 85F 
L5 18L 
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28 


42 30L 






42 31L j 


29 


L5 15L 
42 32L 


f 


30 


46 32L 
50 ( )F 


i 
I 
i 
1 
j 


31 


TJ 7F 






40 ( )F j 


32 






40 ( )F j 


33 


F5 30L 1 




42 30L 


1 


3^ 


42 31L 
L5 32L- 




35 


L4 ii4f 
40 32L 




36 


F5 3F 
40 3F 




37 


LO 17L 
32 30L 




38 


22 38L 
40 3F 




39 


L5 18L 
46 46L 




40 


L5 16l 

42 47L 




4i 


46 kjL 
22 43L 




42 


L5 114F 
L4 46l 




43 


46 46l 
L5 18L 




14 


42 46l 
F5 3F 
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45 


40 3F 

41 2F 


i 




46 


50 ( )F 
7J ( )F 






47 


L4 ( )F 
40 ( )F 






48 


F5 46l 
42 46L 






49 


L4 47L 






50 


40 47L 
F5 2F 






51 


40 2F 




. 


52 


36 46l 






53 


LO 17L 
36 42L 






54 


f4 8f 








42 8F 






55 


LO 97 








32 25L 




Test for last column of partitioned uatrix 


56 


92 63F 








26 111F 




Go to prist out C 




00 671 






B 


outia* f 2 - 52 


Print (A) with ©r viiiwwEt Sign, to a Places 




00 85K 1 


■ 


B 


outine IT 2 * 88 


Input a Sequence of Dfteiail Fractions 




00 1UK 









HI 8F 
4l 1023F 






1 


00 F 








00 10F 
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LOCATION 


OKDEK 


.2 


80 F 




00 4F 


3 


00 IF 




j 00 IF 




j 24 Iks 




j 00 111K 





| hi 31F 




50 5F 


1 


75 5F 




26 10F 




00 10K 





L5 5F 




Sk F 


1 


10 IF 




36 2L 


2 


iA 26F 




k2 3L 


3 


4l 50F 




L5 ( )F 


if 


50 3k 


i 

i 


50 kL 


! 5 


26 67F 


i 
i 


F5 30F 


1 6 


ho 31F 




LO 5F 


7 


36 9L 




F4 3L 


8 . 


42 3L 




22 3L 


9 


92 770F 


i 


92 161F 


; 10 


F5 31F 


■ 


ho 31F 
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Constants 



Stop to read in parameters 
Beginning of Part II of Program tape 



— Print the product matrix by columns 



Punch N 
Punch l.f.'s 



LOCATION 



ORDER 



11 
12 

13 
Ik 



LO 120F 
36 llKL 
F5 31F 
iA 3L 
42 3L 
26 3L 
92 191P 
OF F 
00 TK 
00 F 

00 3162 2776 6500 J 

2k 115N 
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Punch l.f.'s 
End 

— Scaling constant 

Go to clear numerical storage 



